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We present an analysis of the electric dipole moment (EDM) of the deuteron as induced by CP- 
violating operators of dimension 4, 5 and 6 including 6qcd, the EDMs and color EDMs of quarks, 
four-quark interactions and the Weinberg operator. We demonstrate that the precision goal of the 
EDM Collaboration's proposal to search for the deuteron EDM, (1 — 3) x 10~^'^ecm, will provide an 
improvement in sensitivity to these sources of one-two orders of magnitude relative to the existing 
bounds. We consider in detail the level to which CP-odd phases can be probed within the MSSM. 



The most stringent constraints on flavor-diagonal CP 
violation in the hadronic sector arise from bounds on the 
EDMs of the neutron [1], mercury [2], and in certain 
cases thallium [3]. These experiments have important 
implications for physics beyond the Standard Model, and 
its supersymmetric extensions in particular (see e.g. [4]). 

In what follows, we will show that a proposed measure- 
ment of the deuteron EDM [5], with projected sensitivity 

Idcl < (1 -3) X 10"27g(,j^^ 

would improve the sensitivity to Oqcb and SUSY CP- 
violating phases by one to two orders of magnitude. We 
find that the dependence of do on the underlying QCD- 
sector CP-odd sources is closest to dug and is complemen- 
tary to dn- Moreover, in addition to the improvement in 
precision, djj has a significant advantage over dug due to 
the rather transparent nuclear physics in the former and 
thus smaller theoretical uncertainties. Consequently, the 
experiment will be able to probe classes of supersymmet- 
ric models which escape the current EDM bounds. 

We now proceed to analyze the deuteron EDM do, 
defined via the interaction of the deuteron spin / with an 
electric field, Ti. = —do I ■ E, working upwards in energy 
scale. Starting at the nuclear level, the deuteron EDM 
receives contributions from a singlet combination of the 
constituent proton and neutron EDMs, but also arises 
due to meson (predominantly pion) exchange between 
the nucleons with CP-odd couplings at one of the meson- 
nucleon vertices. Thus, we can represent the EDM as 

do = (d„ + dp) + djl"^, (2) 

where the third term includes the meson-exchange con- 
tribution and depends on the CP-odd pion nucleon cou- 
plings, 

^(7P - S^^nnNt^'Ntt'^ + gi'lj^NNir" . (3) 



In a recent analysis, Khriplovich and Korkin [6] (see also 
[7]) showed that dj^'^ receives a dominant contribution 
from the isospin-triplet coupling g^^^ . In a zero-radius 
approximation for the deuteron wavefunction, the result 

^ttNN _ (^g^^NNg^T^NN 1 +^ ^,^ 



depends on the parameter ^ = ^mpe/mT^, determined by 
the deuteron binding energy e — 2.23 MeV. Numerically, 
this implies 

^ -(1.3 ± 0.3) e [GeY-\ (5) 

a result that can be improved systematically, and the 
error correspondingly reduced [6], with the use of more 
realistic deuteron wave functions. 

To make direct contact with models of CP violation, 
we require the dependence of dn, dp, and g'^^^ on the pa- 
rameters in the underlying CP-odd Lagrangian at 1 GeV. 
Up to dimension five, the relevant hadronic operators are 
the ^-term and the EDMs and color EDMs (CEDMs) of 
quarks 

q—u,d,s 

(6) 

where GG = e^,^paGf""'GP'^'' /2 and Ga = t<'G^""'(J^,^. 
Note that the dimension-six Weinberg operator, GGG, 
as well as numerous four-quark operators, may, in certain 
models, also contribute at a similar level to the quark 
EDMs and CEDMs. 

Models of new CP- violating physics can be cast into 
two main categories: (i) models that have no Peccei- 
Quinn (PQ) symmetry [8] and exact CP or P symmetries 
at high energies and consequently = at tree level; and 
(ii) models that invoke a Peccei-Quinn symmetry to re- 
move any dependence of the observables on 6. In models 
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FIG. 1: Contributions to (d), with isospin violation through 

r) — -K mixing. 



of the first type, generated by radiative corrections is 
likely to be the main source of EDMs. 

To determine doid), one may first try to make 
use of the chiral techniques [9] that determine the 9- 
induced pion-nuclcon coupling constant, 5^^r7v(^) = 
m^9f~^{N\uu~dd\N) (where to* — mumd/ {mu + md)), 
and the one loop 0(to^ log TOt) contribution to c?„. It 
is easy to see, however, that dniO) is incalculable within 
this approach because the chiral logarithms exactly can- 
cel in the dn + dp combination, and g^'^^O) = unless 
isospin violating corrections are taken into account. 



The cancellation between dn{d) and dp{9) does not 
hold in general. To calculate doiO) we use leading or- 
der QCD sum-rule estimates which imply [10], 



dn{e) + dp{9) = 

-(2±0.8)7r2(- 



rriN 



GeV/ (IGeV)' 



:m^X(i9, (7) 



where {qcp^i, q)F = eqxFf,u{qq) defines the magnetic sus- 
ceptibility X ~ ~(6 — 9) GeV^^ [11] of the vacuum, re- 
cently computed to be at the upper end of this range, 
X = -A^c/(47r2/2), by Vainshtein [12]. The subleading 
corrections to the sum rule were computed and are of or- 
der 10-15% [10], while the uncertainty in x and freedom 
in the choice of nucleon interpolating current lead to a 
larger overall uncertainty of 30-40% [10] . 

It turns out that despite an additional suppression fac- 
tor, the corresponding contribution to g^^^{9) is not neg- 
ligible and contributes to d^ at approximately the same 
level as (7). To take it into account, we note that isospin 
violation arises predominantly through 77 — vr mixing as 
shown in Fig. 1(a). The inverted diagram of Fig. 1(b) 
provides at most a 10% correction, due primarily to the 
small size of gjjNN and {N\uu — dd\N) relative to g-^NN 
and {N\uu-\-dd—2ss\N). Fig. 1(a) leads to the following 
result: 



rud 



4to„ 



'-{N\uu + dd-2ss\N). (8) 



Combining (7) and (8), we obtain 



do{9) = -eO 



2ti 



2 Xm* (qq) 



(1 GeV) 3 
TO* (mrf - TO„) 



4/. 



{N\uu + dd - 2ss\N) 



(9) 



which numerically takes the form 

dDi9) ~ -e [(3.5 ±1.4)-)- (1.4 ±0.4)] 
xlO-^9 [GeV^^], 



(10) 



using standard quark mass ratios [13], and quark conden- 
sates over the nucleon (see e.g. [14]). The second term 
in (10) arises from the CP-odd pion-nucleon interaction. 

This result is interesting for several reasons. Firstly, if 
the projected experimental sensitivity (1) is achieved, a 
null result for dn will imply 



|6l| < 3 X 10" 



(11) 



which represents an improvement of over an order of mag- 
nitude relative to the best current bound arising from the 
limit on the neutron EDM. We note that the recent in- 
clusion of many-body effects in the nuclear component of 
the calculation of dng [15] has led to a significant reduc- 
tion of d-iig{g^^'>), thus relaxing the mercury EDM con- 
straint on 9 by an order of magnitude. It is also impor- 
tant to note that the two sources for 9 in (9) have quite 
different origins, and thus a cancellation would be unnat- 
ural. Given the relatively good theoretical control over 
the contribution entering through g^-^\ the uncertainty 
in the estimate (7) is of less concern. The bound (11) 
has important implications for solutions to the strong 
CP problem within supersymmetry. In particular, the 
left-right symmetric SUSY models typically predict 9 in 
the range 10~^ — 10"^" [16], allowing a direct probe via 
the d]j experiment. 

Introducing a PQ symmetry allows the axion to re- 
lax to its minimum thereby rendering 9 unobservable. 
Adopting this approach, we are left with the dimension 
five quark EDMs and CEDMs as the leading candidates 
for the position of dominant CP-odd source. The con- 
stituent EDMs of the proton and neutron receive con- 
tributions from both of these operators, with the QCD 
sum-rules result (omitting for now the Weinberg opera- 
tor) [17] 

dn{dq,dq) + dp{dq,dq) ~ (0.5 ± 0.3) ((i„ + drf) 

-(0.6 ± 0.3)e [(J„ - dd) + 0.3(J„ + dd)\ , (12) 

where we have split the CEDM contribution into singlet 
and triplet combinations. A possible contribution from 
ds is removed at this order under PQ relaxation. The 
quoted errors have the same origin as those in (7) for the 
dependence of c?„ and dp on 9. 
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The triplet pion nucleon coupling g^^'> receives a domi- 
nant contribution from the triplet combination {du — dd) 
of CEDMs, and the "best" value for this coupling was 
recently determined using sum- rules [18], 



y-nNN 



5+4 



X 10 



-12 



du — dd 
10-26 cr 



(13) 



with a rather large (overall) uncertainty due to an ex- 
act cancellation at the level of vacuum factorization. We 
quote the non-Gaussian errors determined via parameter 
variation [18]. Since this result enters without any addi- 
tional isospin- violating suppression factor, it numerically 
dominates the CEDM contribution to dn. Combining 
(12) and (13), we find 

doidq, dq) ~ -e{du - dd) [5iji + (0.6 ± 0.3)] 
-(0.2 ± 0.1)e(d„ + dd) + (0.5 ± 0.3)(d„ + dd), (14) 

where the constituent nucleon EDMs provide a 10% cor- 
rection to the triplet CEDM contribution. We con- 
clude from this result that for models with edi ~ di the 
deuteron EDM is predominantly sensitive to the triplet 
combination of CEDMs, as is the mercury EDM. More- 
over, if the predicted precision is achieved, its sensitiv- 
ity to the triplet CEDM combination at the level of a 
fewxlO~28 ecm would represent an improvement on the 
current mercury EDM bound by two orders of magni- 
tude. 

We now turn to an analysis of the predicted sensi- 
tivity to new CP-odd sources focusing on the minimal 
supersymmetric standard model (MSSM) with universal 
boundary conditions at the GUT scale for all parameters 
except for those in the Higgs sector. This exception al- 
lows us to satisfy all phenomenological and cosmological 
constraints for a wide range of squark masses while keep- 
ing the other parameters fixed [19]. In this case, there 
are two CP violating phases, identified with the phases 
of the n parameter in the superpotential and the phase 
of a common trilinear soft-breaking term Aq. 

In Fig. 2, we plot the EDMs as a function of the left- 
handed down squark mass by varying mo from 0.25 - 
10 TeV, while keeping mi/2 (as well as the other input 
parameters) fixed. For this choice of parameters, the light 
Higgs mass is about 120 GeV and the lightest neutralino 
is a mixed gaugino/Higgsino state. The curves begin at 
rhdj^ ~ 1.2 TeV corresponding to mp = 250 GeV with 
mi/2 = 600 GeV. In this figure, the theoretical average 
values of the neutron, thallium and mercury EDMs are 
normalized to their current experimental limits, while d^ 
is normalized to 3 x lO"^^ e cm. The theoretical error 
bands are generally very narrow on these log-scale plots 
and are not shown. For low Oi^(a): the EDMs scale with 
6 and therefore the results for other (small) choices of 
(^fi(A) can be deduced from this figure. We immediately 
see that the projected sensitivity of dn to squark masses 
extends beyond 10 TeV, well beyond that of the existing 
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FIG. 2: The EDMs of the deuteron (black), mercury (green), 
the neutron (red), and thaUium (blue) as a function of the SUSY 
soft breaking scalar mass mo, displayed in terms of the left- 
handed down squark mass. In a) By^ = 0, (9^ = tt/IO and in b) 
8 A = i"/10, 6^ = 0. The EDM is normalized to the experimental 
constraint in each case. 



bounds or the reach of colliders in the foreseeable future. 
Note that the dips observable in the plot of dug for 9f_i ^ 
are due primarily to cancellations between quark CEDM 
and electron EDM contributions. 

The do experiment will also be able to probe a popular 
solution to the SUSY CP problem, the "decouphng" sce- 
nario. This framework assumes that the sfermions of the 
first two generations have masses in the multi-TeV range 
thus suppressing the one-loop EDM contributions to an 
acceptable level and allowing CP-odd phases to be of or- 
der one [20]. To satisfy the cosmological constraints on 
dark matter abundance [21], and to avoid excessive fine- 
tuning in the Higgs sector, the masses of the third genera- 
tion sfermions should be near the electroweak scale. The 
Weinberg operator is then generated at two-loop order, 
providing the primary contribution to d^ [22, 23]: 



dn{w) +dp{w) ^ ew x20 MeV, 



(15) 



where w is the coefficient of the Weinberg operator eval- 
uated at 1 GeV. The Weinberg operator provides a neg- 
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120 GeV (with the same color- 



coding as in Fig. 2). The width of the deuteron band normalized 
to 3 X 10~^^ e cm is too small to be visible on the plot and is 
artificially widened by a factor of 10. 



due to additional chiral 
(1) 



ligible contribution to 

suppression and isospin violating factors in g'^'^pfj^iw) 
Presently, order one CP-violating phases in this frame- 
work are barely compatible with the experimental con- 
straint on dn [24]. Therefore, an improvement in the 
experimental precision by a factor of 10 or more, to the 
level of 10~^^ e cm, would provide a crucial test for these 
models. Failure to observe do would necessarily imply 
that the CP-violating phases are small contrary to the 
primary assumptions of the model. 

Next, we analyze constraints on the SUSY CP- 
violating phases 9a, with the superpartner mass scales 
fixed as shown in Fig. 3. This is a CMSSM point (the 
Higgs soft masses are unified with other sfermion masses) 
with a relatively low Higgs mass of 114 GeV [19]. We 
observe that du combined with the thallium constraint 
can put tight bounds on both phases including 6a that 
is otherwise poorly constrained. An improvement of the 
bound on the triplet CEDM combination by a factor of 30 
or more would allow one to probe SUSY CP-odd phases 
of size 10^^ or below (10~^ or so for the A-terms). In 
a number of theoretically motivated scenarios, phases of 
this size are naturally expected. In particular, if the A- 
terms are hermitian at the GUT scale as happens in the 
left-right and other models, RG running induces small 
phases in the diagonal elements. For a variety of tex- 
tures, the CEDMs of the light quarks are of order 10~^^ 
cm [25], and thus observable at the upcoming experiment. 

Finally, we consider the sensitivity of do to the di- 
mension 6 operators, CijqiQiqjij^qj, which may be im- 
portant in two Higgs doublet models, left-right symmet- 



ric models, and certain supersymmetric scenarios. Typi- 
cally, Cij can be parametrized as Cij = cY^^Y^^^M^^, 
where Y^^ are the SM quark Yukawa couplings, Mh is 
the mass of the (lightest) Higgs boson, and the coefficient 
c is model dependent. Existing EDM bounds are sensi- 
tive to Cij only with the help of an enhancement at large 
tan/3, c ~ tan^ f3 or tan^ /3 [26], or in the top quark sec- 
tor where Ctq induces w and/or light quark (C)EDMs via 
the Barr-Zee mechanism [27]. The projected sensitivity 
to di) would in contrast probe C'ij for all quark flavors 
down to c ~ 0.01 — 0.1 for Mh ~ 100 GeV, thus providing 
valuable constraints even for tan/3 ~ 0(1). 

In conclusion, we have presented an analysis of the 
deuteron EDM in terms of the relevant Wilson coeffi- 
cients and studied the implications of a djj measurement 



at the level of a fewxlO 



;cm. We have shown that 



this would lead to a factor of 10 to 100 gain in sensitivity 
to various CP violating sources of dimension 4, 5 and 6. 
This has important consequences for supersymmetry and 
other scenarios for physics beyond the Standard Model. 
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